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We study a class of discrete-time quantum walks with inhomogeneous coins defined in [Y. Shikano 
and H. Katsura, Phys. Rev. E 82, 031122 (2010)]. We establish symmetry properties of the 
spectrum of the evolution operator, which resembles the Hofstadter butterfly. 

PACS numbers: 03.65.-w, 71.23.An, 02.90.+p 

o: 

Throughout this paper, we focus on a one-dimensional discrete time quantum walk (DTQW) with two-dimensional 
coins. The DTQW is defined as a quantum-mechanical analogue of the classical random walk. The Hilbert space of 
the system is a tensor product Hp® He, where Tip is the position space of a quantum walker spanned by the complete 
orthonormal basis \n) (n 6 Z) and H c is the coin Hilbert space spanned by the two orthonormal states \L) = (1,0) T 
and \R) = (0, 1) T . Here, the superscript T denotes matrix transpose. A one-step dynamics is described by a unitary 
CN ' operator U = WC with 
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C = J2[( a n\n,L} + Cn\n,R}){n,L\ + (d n \n,R} + bn\n,L}){n,R\}, (1) 

n 

W = Y^ i\n - 1, L){n, L\ + \n + 1, R)(n, R\) , (2) 



where |n, £) =: |n) ® |£) £ ri p ® He (£ = L,R) and the coefficients at each position satisfy the following relations: 
\a„.\ 2 + \c„.\ 2 = 1, a n b n +c n d n = 0, c„ = -A n b n , d n = A„a„, where A n = a n d n - b n c n with |A„| = 1. Two operators 
C and W are called coin and shift operators, respectively. The probability distribution at the position n at the tth 
step is then defined by 

Pr(n;t)= ]T \(n,Z\U%ct>)\ 2 . (3) 

■ A homogeneous version of this DTQW was first introduced in Ref. 
Suppose that the coin operator is given by 

C(a, 0) =y^ [{cos{2nan + 2n9)\n, L) + sin(27rcm + 2tt9)\ji, R))(n, L\ 



+(cos(27ran + 2n0)\n, R) - sin(27rcm + 2w6)\n, L))(n, R\] 
:=^2\n){n\»d n {a,e), (4) 

n 

where a and 9 are constant real numbers. Then this class of DTQW is called an inhomogeneous quantum walk (QW). 
This model is based on the idea of the Aubry- Andre model [2J , which provides a solvable example of metal-insulator 
transition in a one-dimensional incommensurate system. In this class of DTQW, we have obtained the weak limit 
theorem as follows. 

Theorem 1 (Shikano and Katsura Q). Fix 9 = 0. For any irrational a £ M\Q and any special rational a=^6Q 
with relatively prime P (odd integer) and Q, the limit distribution of the inhomogeneous QW is given by 

(t->oo), (5) 
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FIG. 1: Probability distribution of the inhomogeneous QW at 300th step with a = 1/3. From simple algebra, it can be easily 
shown that the inhomogeneous QW is finitely confined when 9 — (2m — 1)/12 (m G Z). 



where X t is the random variable for the position at the t step, "=$■" means the weak convergence, and r/ (> 0) is an 
arbitrary positive parameter. Here, the limit distribution I has the probability density function f{x) — 5{x) (x G K), 
where S(-) is the Dirac delta function. This is called a localization for the inhomogeneous QW. 

However, in the case of the other rational a, it has not yet been clarified whether the inhomogeneous QW is localized 
or not. This is still an open question. The situation becomes more complicated when we consider a nonzero 9. As 
seen in Figure [TJ the reflection points for the quantum walker (see more details in Ref. |3j, Lemma 1 and Figure 2]) 
are changed by the parameter 9. In the rest of the paper, we will establish symmetry properties of the eigenvalue 
distribution of the one-step evolution operator (U — WC) at 9 = 0. 

Theorem 2. For the eigenvalues of the one-step evolution operator WC, the following properties hold: 

(PI) All the eigenvalues at a are identical to those at 1 — a. 

(P2) For every eigenvalue X, there is an eigenvalue A*. 

(PS) For every eigenvalue X, there is an eigenvalue —A. 

(P4) All the eigenvalues are simple, i.e., nondegenerate. 

(P5) There are four eigenvalues X = ±1, ±i for any a = G Q. 

(P6) Every eigenvalue X at a = G Q corresponds to an eigenvalue iX at a + 1/2. 

Proof. The proofs of properties (PI) - (P5) can be found in Ref. Q. Here, we give a proof of (P6). According to 
Ref. [3, Theorem 3], the eigenvalues of WC and WC are identical. Therefore, we only study the eigenvalues of CW. 
First, we can express the wavefunction at the tth step evolving from the state |0, 0) by CW: 






(6) 



The one-step time evolution of the coefficients ipt(n, £) is given by 




(7) 



Here, we define t by ipt(n, £) and a square matrix of order 4Q, denoted as CW, as 



t +i = CWi^f, 



(8) 
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FIG. 2: Eigenvalue distribution of the one-step operator for the inhomogeneous QW (U). Arguments of the eigenvalues of WC 
(vertical axis) are plotted as a function of the parameter a — ^ (horizontal axis) with Q < 60. Here, P (odd number) and Q 
are relatively prime. 



see more details in Ref. [3j. Let — (<p(—Q; R), <p(—Q + 1; L), <p(—Q + 1; R), <f(Q', L)) T be the eigenvector of CW 
at a with the eigenvalue A and (p — (ip(—Q;R),cp(—Q + l;L),<p(—Q + 1; R), <p(Q; L)) T be at a + 1/2 with the 
eigenvalue A. Then, according to Eq. (O, we obtain 

\<p(-Q;R) = (-l)^<p(-Q + l;L), 
A^(g;L) = (-l) £ f i ^(Q-l;i?) (9) 

and 

m-Q;R) = (-l)- Q (-l)^^(-g + l;L), 

A^(Q;L) = (-l) c ?(-l)^ i ^(Q-l;i?), (10) 

where we have used the fact C n (a + 1/2) = (— l) n C n (a, 0). Now we apply the following local unitary transformation 
to Eq. (HOI): 

jV(n;£) when n is even, 
I (n;4) when n is odd. 

According to Eq. ©, defined by Eq. ([TTj) can be taken as the eigenvector CW at a with the eigenvalue A = iX. □ 

Figure [2] shows the numerically obtained spectrum of CW as a function of a, which is quite similar to the Hofstadter 
butterfly [J|. By combining all the properties of (P1)-(P6), the smallest fundamental domain of this diagram is 
identified as the triangular region shown in Figure [2] Therefore, we have rigorously established all the symmetries in 
Figure [2] 
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